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ABSTRACT 

We present new HST ultraviolet color-magnitude diagrams of 5 massive Galactic globular clusters: 
NGC 2419, NGC 6273, NGC 6715, NGC 6388, and NGC 6441. These observations were obtained to in- 
vestigate the "blue hook" phenomenon previously observed in UV images of the globular clusters oj Cen and 
NGC 2808. Blue hook stars are a class of hot (approximately 35,000 K) subluminous horizontal branch stars 
that occupy a region of the HR diagram that is unexplained by canonical stellar evolution theory. By coupling 
new stellar evolution models to appropriate non-LTE synthetic spectra, we investigate various theoretical ex- 
planations for these stars. Specifically, we compare our photometry to canonical models at standard cluster 
abundances, canonical models with enhanced helium (consistent with cluster self-enrichment at early times), 
and flash-mixed models formed via a late helium-core flash on the white dwarf cooling curve. We find that 
flash-mixed models are required to explain the faint luminosity of the blue hook stars, although neither the 
canonical models nor the flash-mixed models can explain the range of color observed in such stars, especially 
those in the most metal-rich clusters. Aside from the variation in the color range, no clear trends emerge in the 
morphology of the blue hook population with respect to metallicity. 

Subject headings: globular clusters: general - globular clusters: individual (NGC 2419, NGC 6273, NGC 6715, 
NGC 2808, NGC 6388, NGC 6441) - stars: atmospheres - stars: evolution - stars: hori- 
zontal branch - ultraviolet: stars 



1. INTRODUCTION 

Understanding the formation history and chemical evo- 
lution of the Galactic globular clusters has always been a 
formidable, but fundamental, task. However, this task has 
become even more challenging with the recent discovery of 
complex stellar populations within individual globular clus- 
ters. The most well-known example is cj Cen (Anderson 
1997; Bedin et al. 2004; Ferraro et al. 2004), but the same 
phenomenon has also been discovered in other massive glob- 
ular clusters, including NGC 2808 (D' Antona et al. 2005; Pi- 
otto et al. 2007), NGC 1851 (Milone et al. 2008; Han et al. 
2009), NGC 6715 (M54; Layden & Sarajedini 2000; Siegel 
et al. 2007), and NGC 104 (47 Tuc; Anderson et al. 2009). 
In Lj Cen, accurate photometry demonstrates that its main se- 
quence (MS) population is split into multiple sequences, with 
the bluest MS stars being more metal rich, implying that the 
corresponding sub-population is perhaps enhanced in helium 
(Y ^ 0.4; Piotto et al. 2005). Such helium-rich stars may have 
formed in a second stellar generation from the helium-rich 
ejecta of the initial stellar generation that formed the clus- 
ter (e.g., see Renzini 2008). Possible candidates for produc- 
ing these ejecta include asymptotic-giant-branch (AGB) stars, 
massive MS stars, type II supernovae, and even stars external 
to the cluster (Bekki & Norris 2006 and references therein). 
These clusters thus refute the long-standing premise that glob- 
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ular clusters are simple stellar populations. 

The discovery of complex populations in globular clusters 
has understandably generated much excitement, especially in 
regard to the morphology of the horizontal branch (HB). Be- 
cause the MS turnoff (MSTO) mass decreases strongly with 
increasing helium at a fixed cluster age, a high helium abun- 
dance will lead to a bluer HB morphology for a given range 
in the red-giant branch (RGB) mass loss, and this, in turn, 
might account for the long blue HB tails found in some clus- 
ters (D'Antona et al. 2002; Busso et al. 2007). In this sce- 
nario, the extreme HB (EHB) stars, located at the hot end of 
the HB at temperatures Tgff > 16,000 K and surface gravities 
log ,g > 5, would be the progeny of the most helium-rich MS 
stars. The hottest EHB stars have extremely thin envelope 
masses (^ lO"-' Mq), due to extensive mass loss on the RGB. 
The analogs of the EHB stars in the field are the subdwarf 
B (sdB) stars. Both unresolved spectroscopy and resolved 
imaging have demonstrated that EHB stars are the dominant 
contributors to the "UV upturn" in the otherwise cool spectra 
of elliptical galaxies (Brown et al. 1997; Brown et al. 2000; 
Brown et al. 2008). 

UV observations of globular clusters have revealed yet an- 
other puzzling result, namely, that two of the massive clus- 
ters with evidence for a helium-rich MS population (lj Cen 
and NGC 2808) also contain an unexpected population of 
subluminous EHB stars that cannot be explained by canon- 
ical stellar evolution theory. These subluminous stars were 
first discovered in oj Cen (D'Cruz et al. 1996, 2000), where 
they form a "blue hook" (BH) feature at the hot end of the 
canonical EHB. D'Cruz et al. (1996) proposed that such stars 
may have undergone a delayed helium-core flash after peeling 
away from the RGB due to high mass loss. 

We now know that stars can follow a variety of evolutionary 
paths to the HB, depending on when the helium-core flash oc- 
curs (see Figure 1). For normal rates of mass loss, the helium- 
core flash will occur at the tip of the RGB (RGB tip flasher. 
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Figure 1. Various evolutionary paths for producing an HB star. The zero-age HB (ZAHB) phase is highlighted in grey, while the pre-ZAHB evolution {solid 
curves) and post-ZAHB evolution {dashed curves) are in black. The peak of the helium-core flash is marked by an asterisk. The four panels show the evolution for 
progressively larger amounts of mass loss on the RGB. The evolution in the first two panels produces canonical HB stars in which the H-rich surface composition 
does not change during the helium-core flash. In the third panel, the helium-core flash occurs on the WD cooling curve, producing a flash-mixed star having 
a surface composition highly enriched in helium and carbon and a temperature significantly hotter than the canonical HB. In the fourth panel, the helium flash 
never occurs and the star dies as a helium WD. 



Figure la). However, for sufficiently high rates of mass loss, 
the helium-core flash can occur either during the crossing of 
the HR diagram (early hot flasher. Figure lb) or on the white 
dwarf (WD) cooling curve (late hot flasher, Figure Ic), as first 
shown by Castellani & CasteUani (1993). Brown et al. (2001; 
hereafter Paper I) demonstrated that a delayed helium-core 
flash on the WD cooling curve would result in a stellar atmo- 
sphere extremely enhanced in helium (^96% by mass) and 
carbon (~4% by mass), due to the mixing of the stellar en- 
velope into the hot heUum-buming core. This "flash mixing" 
decreases the opacity below the Lyman limit, thereby lower- 
ing the flux at longer wavelengths, and increases the effec- 
tive temperature in the EHB stars. In Paper I, we found that 
both of these effects together might explain the BH feature in 
the color-magnitude diagram (CMD) of NGC 2808. Evidence 
for BH populations in other clusters soon followed, including 
NGC 6715 (Rosenberg et al. 2004), NGC 6388 (Busso et al. 
2007), NGC 2419 (Ripepi et al. 2007), and NGC 6273 (noted 
by Rosenberg et al. 2004, but originally appearing in Piotto et 
al. 1999). NGC 6441 may host BH stars, although the case 
has been unclear (Busso et al. 2007; Dieball et al. 2009). 

To extend the investigations enabled by UV observations of 
ij Cen and NGC 2808, we have obtained far-UV and near- 
UV images of those additional clusters with the strongest ev- 
idence for BH populations: NGC 2419, NGC 6273 (M19), 
NGC 6715 (M54), NGC 6388, and NGC 6441. A charac- 
terization of BH populations in globular clusters can only be 
performed definitively with an ultraviolet CMD, where, due 
to the small bolometric correction, the photometric effects of 
enhanced helium and carbon abundances can be distinguished 
from changes in temperature. These clusters span a wide 
range of metallicity (-2.1 < [Fe/H] < -0.5), which enables 
an exploration of how the BH population is distributed in tem- 
perature, luminosity, and size, relative to the canonical EHB 
population, as a function of cluster environment and metaUic- 



ity. 

To avoid confusion, it is worth clarifying our terminology. 
BH stars were first defined observationally as being hotter and 
fainter than the hot end of the canonical EHB. However, we 
will show in this paper that some clusters host stars fainter 
than the canonical EHB, but with somewhat redder colors. 
These red colors may imply temperatures down to ^15,000K. 
The stars appear distinct from the still redder blue straggler 
sequence (BSS), and they appear to be associated with the 
bluer BH stars. We will also consider these to be BH stars, 
without necessarily implying that they have the same physical 
origin. 

In this paper, we compare our UV CMDs of five mas- 
sive globular clusters, together with the previous UV CMD 
of NGC 2808, to the expectations from canonical evolution 
theory and flash-mixing. In each case, we explore the ex- 
pectations for stars bom at standard cluster abundances and 
for stars born with enhanced helium abundances. To ensure 
a consistent exploration of these scenarios, our evolutionary 
models are tailored to the appropriate abundance profiles, as 
are the model atmospheres and synthetic spectra employed to 
transform from the theoretical plane to the observable plane. 

2. OBSERVATIONS AND DATA REDUCTION 

The six clusters discussed herein have all been imaged in 
the far-UV and near-UV using the Hubble Space Telescope 
(HST). Due to variations in the instrumentation available for 
each observation, the bandpasses were not identical for each 
cluster (see Figure 2), but the effective wavelengths in each 
wavelength region (far-UV or near-UV) are similar, easing 
comparison of CMDs constructed from these bandpasses. In 
Table 1, we summarize the observations for each cluster We 
also include NGC 2808, the subject of Paper I, observed with 
the Space Telescope Imaging Spectrograph (STIS); each STIS 
UV channel has pixels of 0.025 x 0.025 arcsec and a field of 
view of 25 x 25 arcsec. 
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Table 1 

Ultraviolet Imaging 



Cluster 


t'ar-UV 




near-U V 






bandpass 


exposure (s) 


bandpass 


exposure (s) 


NGC2419 


ACS/SBC/F150LP 


5180 


ACS/HRC/F250W 


5080 


NGC6273 


ACS/SBC/F150LP 


1980 


WFPC2A¥F3/F255W 


11000 


NGC6715 


ACS/SBC/F150LP 


1980 


WFPC2A¥F3/F255W 


11000 


NGC2808 


STIS/FUV/F25QTZ 


8916 


STIS/NUV/F25CN270 


8906 


NGC6388 


ACS/SBC/F150LP 


1980 


WFPC2A¥F3/F255W 


11000 


NGC6441 


ACS/SBC/F150LP 


4680 


WFPC2A¥F3/F255W 


11000 



As originally proposed, the observations in the current pro- 
gram intended to use the Solar Blind Channel (SBC) and High 
Resolution Channel (HRC) of the Advanced Camera for Sur- 
veys (ACS) on HST, with the SBC F150LP filter providing 
the far-UV bandpass and the HRC F250W filter providing the 
near-UV bandpass. The ACS/SBC has pixels of 0.034 x 0.030 
arcsec and a field of view (FOV) of 34.6 x 30.1 arcsec, while 
the ACS/HRC has pixels of 0.028 x 0.025 arcsec and a FOV 
of 29 X 26 arcsec. Due to its larger distance and smaller angu- 
lar size, the observations of NGC 2419 were planned with a 
single tile imaged with the SBC and HRC, while the remain- 
ing four clusters were planned with a 2 x 2 mosaic of tiles in 
each bandpass. Due to guide star problems and instrument 
failures, the program did not proceed according to plan. 

In December 2006, the first observations of NGC 2419 were 
attempted. The HRC exposures were obtained in two con- 
secutive orbits, with no problems. The SBC exposures were 
obtained in the subsequent two orbits, but with serious prob- 
lems. The guide star re-acquisition failed in the first SBC 
orbit, setting the Take Data Flag down and resulting in two 
blank images. In the final orbit, the re-acquisition of guide 
stars was successful, but due to a flight software bug, the tele- 
scope did not complete the maneuvers needed to point at the 
correct field. As a consequence, the next two exposures were 
not blank but offset by nearly an arcmin from the success- 
ful HRC images. This flight software bug, which was ap- 
parently unique to the HST two-gyro mode in the situation 
where successful guide star acquisitions followed guide star 
failures in the same visit, has since been fixed. The SBC 
exposures for NGC 2419 were eventually obtained a year 
later. In short, NGC 2419 was observed with the intended 
ACS/SBC/F150LP and ACS/HRC/F250W bandpasses, after 
some difficulties. 

Due to a failure of an ACS power supply in January 2007, 
the HRC channel stopped working. At that time, only the 
near-UV (HRC) images of NGC 2419 had been taken, so the 
near-UV images of the remaining four clusters were obtained 
with the Wide Field and Planetary Camera 2 (WFPC2), us- 
ing the F255W filter and the WF3 chip. The WF3 chip has 
0.1 X 0.1 arcsec pixels and a FOV of 80 x 80 arcsec. The 
near-UV channel on WFPC2 is not as sensitive as that on the 
ACS/HRC, but the larger field of view allowed imaging of the 
intended area in a single pointing instead of the 2x2 mosaic 
of HRC images. Thus, NGC 6273, NGC 6388, NGC 6441, 
and NGC 6715 were observed with the ACS/SBC/F150LP 
and WFPC2/F255W bandpasses. 

Due to these distinctions, the observations of NGC 2419 
and the remaining four clusters were processed somewhat dif- 
ferently. For NGC 2419, the images in each channel were reg- 
istered and drizzled to a common scale of 0."025, while for 
the other four clusters, the scale was 0."0455. The ACS/SBC 
images were combined using the drizzle package (Fruchter 
& Hook 2002), but required no cosmic-ray rejection because 
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Figure 2. The bandpasses employed in our UV imaging (see Table 1) 



the SBC is a photon-counting multianode microchannel array 
(MAMA). The WFPC2 and ACS/HRC detectors are CCDs, 
so images from these detectors were combined with masking 
of cosmic rays. 

Because globular clusters are much sparser in the UV than 
in the optical, we did not perform photometry using fitting 
of the point spread function. Instead, we performed aper- 
ture photometry in each channel with an aperture correction 
to an infinite aperture, using the standard aperture correc- 
tions for these instruments (Heyer et al. 2004; Maybhate et 
al. 2010; Sirianni et al. 2005). For the four clusters observed 
with the ACS/SBC and WFPC2, the apertures were 2.5 pix- 
els and 5 pixels, respectively; for the cluster observed with 
the ACS/SBC and ACS/HRC (NGC 2419), the apertures were 
2.5 pixels and 1.5 pixels, respectively. Because it is a MAMA, 
the SBC requires no correction for charge transfer inefficiency 
(CTl), but the near-UV photometry was corrected for CTl us- 
ing the corrections of Dolphin (2000) and Chiaberge et al. 
(2009), for the WFPC2 and ACS/HRC, respectively. 

The CMD for each cluster is shown in Figure 3. In Figure 
4, we highlight one of these clusters (NGC 6715) to provide 
a guide to the various classes of stars discussed in this pa- 
per: blue HB (BHB) stars, EHB stars, BH stars, BSS, and 
the AGB-Manque (AGBM) progeny of the EHB stars. Note 
that we chose NGC 6715 not because it is particularly repre- 
sentative of massive clusters but because each class of star is 
well-populated in this CMD. Each cluster exhibits a distinct 
hook feature at the blue end of its HB locus. As discussed ex- 
tensively in Paper 1, this hook feature cannot be due to photo- 
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metric scatter, differential reddening, or instrumental effects, 
because these mechanisms could not affect the photometry of 
the EHB without also affecting the photometry of the neigh- 
boring BHB stars, which are only slightly redder; the EHB 
and BHB stars have similar far-UV luminosity and photomet- 
ric errors. Despite the similar photometric errors in the EHB 
and BHB stars, the EHB locus spans a luminosity range that 
is several times larger than that on the hot end of the BHB. 

3. MODELS 
3.1. Stellar Evolutionary Models 

The stellar evolutionary models used in the present study 
were computed with a procedure similar to the one described 
in Paper I. For each cluster, we constructed grids of evolu- 
tionary sequences that followed the evolution of appropriate 
low-mass stars continuously from the zero-age MS through 
the end of the HB phase for two values of the initial helium 
abundance: a helium-normal abundance Y of 0.23, represent- 
ing the helium abundance in stars born during the first stellar 
generation with standard cluster abundances, and a helium- 
rich abundance Y of 0.40, representing a reasonable upper 
limit on the helium enhancement in the second stellar gen- 
eration born after a process of cluster self-enrichment. The 
hotter HB sequences in these grids were also evolved through 
the post-HB phases. 

The heavy element abundance Z for these sequences was 
determined from the observed [Fe/H] value for each clus- 
ter, assuming an a-element enhancement of [a/Fe] = 0.3, 
as appropriate for metal-poor globular cluster stars (Car- 
ney 1996). Using the prescription given by Salaris et al. 
(1993) for converting an a-enhanced abundance into an 
equivalent scaled-solar abundance, we obtained scaled-solar 
heavy element abundances Z of 0.00023, 0.00064, 0.00081, 
0.0015, 0.0077, and 0.0090, respectively, for the six clusters 
NGC 2419, NGC 6273, NGC 6715, NGC 2808, NGC 6388 
and NGC 6441. We used the same Z value to compute both 
the helium-normal and helium-rich sequences for each clus- 
ter. The initial main sequence mass for each composition (Y, 
Z) was calculated for an age of 13 Gyr at the tip of the RGB, 
although ages a few Gyr younger would have little qualitative 
effect on the subsequent evolutionary behavior 

The sequences for each composition {Y, Z) differ from each 
other only in the amount of mass loss along the RGB, which 
we parameterize with the Reimers (1975, 1977) mass-loss for- 
mulation: 

M = -4 X lO-^^r]RL/gR (Mq yr"'), 
where L, g, and R are the stellar luminosity, gravity, and ra- 
dius, respectively, in solar units, and rjR is the well-known 
Reimers mass-loss parameter Because our calculations fol- 
low the evolution through the helium-core flash, we are able 
to determine how the evolutionary path from the MS to the 
ZAHB depends on the assumed mass-loss rate (see Figure 1). 
The red end of the ZAHB is set by the sequence with rju = 0.0 
(no mass loss). Up to a given value of rjR (~0.7 for interme- 
diate [Fe/H] values), the models remain tightly bound to the 
RGB until the helium-core flash (Figure la). However, for 
higher values of the models leave the RGB and evolve to 
high effective temperatures before igniting helium (Castellani 
& Castellani 1993), either during the crossing of the CMD 
(Figure lb) or while descending the WD cooling curve (Fig- 
ure Ic). D'Cruz et al. (1996) suggested that such "hot flash- 
ers" might provide another channel for populating the EHB. 
Although the masses of the late hot flashers for each com- 



position (Y, Z) fall within a narrow range, they are produced 
over a rather wide range in the mass loss rate (Ar/j? ~0. 1). Ex- 
tremely fine tuning of the mass-loss rate is thus not required to 
produce a late hot flasher (see Paper I). For sufficiently high 
values of rju, the models simply cross the HR diagram and 
descend the WD cooling curve without igniting helium, thus 
dying as helium WDs (see Figure Id). 

The RGB tip flashers and early hot flashers ignite helium 
at a time when the models are very luminous. The strong 
hydrogen-burning shell in these models forms a high entropy 
barrier that prevents the flash convection from penetrating 
into the envelope, and hence no mixing then occurs between 
the envelope and the core. Thus any star that undergoes 
the helium-core flash before reaching the WD cooling curve 
will arrive on the ZAHB via canonical evolution (i.e., with 
a normal atmospheric composition). In contrast, stars that 
ignite helium on the WD cooling curve have much weaker 
hydrogen-burning shells and thus much lower entropy barriers 
(Sweigart 1997; Paper I). Under these conditions, penetration 
of the flash convection into the envelope becomes inevitable. 
This flash mixing will mix hydrogen-rich material from the 
envelope into the core and helium- and carbon-rich material 
from the core into the envelope, thereby greatly enhancing 
the surface helium and carbon abundances. A similar pro- 
cess may produce other extremely hydrogen-poor stars (also 
known as "born again" stars) during a very late helium-shell 
flash (Renzini 1990). 

Due to the numerical difficulties associated with this mixing 
process (Paper I), we terminated our sequences as soon as they 
encountered flash mixing. Nevertheless, we can still predict 
the surface composition that these flash-mixed models should 
have when they arrive on the ZAHB. The available calcula- 
tions clearly indicate that flash mixing will consume most, if 
not all, of the envelope hydrogen. Thus we expect the surface 
composition of the flash-mixed models to be strongly depleted 
in hydrogen and enriched in helium. The surface carbon abun- 
dance should also be close to its value in the flash convection 
zone, namely, 4% by mass. For these reasons the flash-mixed 
models given in this paper will assume that the surface carbon 
abundance has been increased to 4% by mass and that all of 
the remaining surface hydrogen has been burned to helium. It 
is worth noting that Cassisi et al. (2003) and Miller Bertolami 
et al. (2008) were able to follow the flash mixing process in 
detail, and their values for the flash-mixed surface composi- 
tion are in good agreement with the composition adopted here. 
Note also that some of the carbon may be burned to nitrogen 
during the flash mixing, depending on the details of the nu- 
cleosynthesis; such an outcome would not have a significant 
effect on the stellar evolution models. 

3.2. Non-LTE Model Atmospheres and Synthetic Spectra 

In order to produce theoretical spectra for the evolutionary 
tracks, we calculated a series of 60 non-LTE line-blanketed 
model atmospheres with our model atmosphere program, 
TlusTY, version 201^. For each globular cluster, we used the 
appropriate metallicities and helium abundances for the vari- 
ous evolutionary scenarios explored. Tlusty computes stel- 
lar model photospheres in a plane-parallel geometry, assum- 
ing radiative and hydrostatic equilibria. Departures from LTE 
are explicitly allowed for a large set of chemical species and 
arbitrarily complex model atoms, using our hybrid Complete 
Linearization/Accelerated Lambda Iteration method (Hubeny 

' Available at |http://nova.astro.umd.edu| 
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Figure 3. Ultraviolet CMDs for six massive globular clusters spanning a wide range in [Fe/H] {points), along with the HB locus expected from canonical stellar 
evolution theory at standai'd helium abundance (Y = 0.23; violet shaded area). The assumed cluster pai'ameters are indicated (grey labels). An approximate 
conversion (assuming normal cluster abundances) between observed color and effective temperature is shown on the upper abscissa in each panel. Photometric 
errors (blue bars) are only significant for stars much redder and fainter than the EHB stars that are the focus of this paper The BH population appears as a 
downward hook at the blue end of the observed HB, where it deviates below the canonical EHB. For the two metal-rich clusters, that deviation occurs far to the 
red of the hot end of the canonical HB. 
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Figure 4. The ultraviolet CMD of NGC 6715 (points), along with the ex- 
pected HB locus based upon canonical evolutionary models at standard he- 
lium abundance (violet shaded area). The various classes of stars are indi- 
cated (colored Imxes). The BH stars are the subluminous subset of the EHB 
stars. Reddening vectors for both redder and bluer stars are indicated (labeled 
grey arrows) 

& Lanz 1995). More specifically, the model atmospheres al- 
low for departures from LTE for 1132 levels and superlevels 
of 52 ions, H I, He I, He II, CI - CIV, N I - N V, O I - O VI, 
Ne I - Ne IV, Mg II, Al II, Al III, Si II - Si IV, P IV, P V, S II 
- S VI, Fe II - Fe VI. Details of the model atom setup are pro- 
vided in Lanz & Hubeny (2003, 2007), and in Cunha et al. 
(2006) for updated Neon models. 

The model atmospheres span the range of parameters suit- 
able for most EHB stars: 20,000 K< Teff <50,000K (with 
a 5,000 K step), 4.5 < log ^ < 6.0 (0.5 dex step), and 
Viurb = 5 km/s. We assumed three different helium composi- 
tions: standard helium (Y = 0.23), helium-rich {Y = 0.4), and 
flash-mixed {Y = 0.96, with carbon and nitrogen mass frac- 
tions of 3% and 1%, respectively; Lanz et al. 2004 - hereafter 
Paper II). Note that the assumption of 4% carbon or 4% ni- 
trogen (or some intermediate mix) makes little difference in 
broad-band UV colors (see Paper 1). 

With our spectrum synthesis code, Synspec, we then cal- 
culated the detailed emergent ultraviolet spectrum of each 
model atmosphere. Quadratic interpolations in the spectra 
database were performed to derive theoretical far-UV and 
near-UV magnitudes for the specific stellar parameters along 
each evolutionary track. The parameters of the quadratic in- 
terpolation are wavelength-dependent. Detailed tests show 
that this interpolation process yields colors that do not differ 
by more than 0.01 mag, compared to using a model atmo- 
sphere calculated with the exact stellar parameters. This ap- 
proach saves significant computing time without significantly 
affecting the final accuracy. 

3.3. Transfer of Models to the Color-Magnitude Diagrams 

With these models in hand, we transformed the evolution- 
ary models for each cluster to the observable plane. For each 
cluster, we interpolated in the grid of synthetic spectra hav- 
ing the appropriate abundance pattern (cluster metallicity with 
Y = 0.23, cluster metallicity with Y = 0.4, or the flash-mixed 
composition). For stars cooler than the lowest effective tem- 
perature in our non-LTE grids, we interpolated in the Castelli 
& Kurucz (2003) spectral grids at enhanced [a/Fe]. The in- 



terpolated spectrum for each stellar evolutionary model was 
reddened using the average Galactic extinction curve of Fitz- 
patrick (1999), and then passed through the relevant near- 
UV and far-UV bandpasses using the IRAF SYNPHOT pack- 
age, resulting in magnitudes in the STMAG system, where 
m = -2.5x logio/A-21.1 mag. Note that this method appro- 
priately applies the reddening in a bandpass-dependent and 
Teff-dependent manner. 

4. COMPARISON OF PHOTOMETRY TO MODELS 
4. 1 . Canonical Models at Standard Helium Abundance 

We first compare our CMDs to the expectations from 
canonical stellar evolutionary models with a standard helium 
abundance. We use the BHB as a standard candle to guide the 
normalization of these models. Specifically, we adopt metal- 
licity and distance values that are representative of values in 
the literature, and then make small adjustments to the extinc- 
tion so as to place the BHB of a canonical model at the ex- 
pected location. Note that the reddening vector in our nifuv 
vs. nifuv-niNuv CMDs is nearly vertical, so one could ad- 
just the distance instead of the reddening, or some combina- 
tion thereof; a change of 0.01 mag in E{B-V) is equivalent 
to a change of 0.079 mag in distance modulus. For the four 
metal-poor clusters, the theoretical HB distribution should be 
coincident with the observed HB locus at the BHB. For the 
two metal-rich clusters. Rich et al. (1997) have shown that 
the HB slopes strongly upward as one moves from the red- 
dest HB stars to the BHB (see also Busso et al. 2007), such 
that the BHB is ^^0.5 mag brighter in the optical than one 
would expect for a canonical locus normalized to the red HB. 
It would be simpler to consistently align our models to the 
red HB in all of our clusters, but it falls outside of our UV 
CMDs, driving our use of the BHB; thus, for the two metal- 
rich clusters we choose an extinction that places the canon- 
ical locus 0.5 mag fainter in the far-UV than the observed 
BHB locus. We are assuming that a star at a given reference 
point on the BHB (defined at a particular effective tempera- 
ture) has changed in bolometric luminosity, such that the off- 
set in the far-UV is the same 0.5 mag observed in the opti- 
cal. The physical motivation for this assumption comes from 
the fact that the bolometric luminosities of all HB stars in 
the range 6000 < Teff < 18000 K are increased by nearly the 
same amount when Y is enhanced, and that the effect on the 
UV-to-optical color in the emergent spectrum is minimal. As 
we shall see in the next section, our choice of extinction is 
equivalent to aligning canonical BHB models of intermediate 
Y (--0.3) to the observed BHB locus. 

In Figure 3, we show the CMD of each cluster compared to 
the expected locus for stars in the core helium-burning stage, 
assuming standard cluster metallicity and helium abundance. 
The comparison is made over a wide color range that includes 
stars much redder than the EHB. In Figure 5, we restrict the 
comparison to those stars hotter than ~9,000 K, and shade 
the HB locus to indicate the timescales expected for HB evo- 
lution. We expect most of the HB stars to follow close to the 
low-luminosity edge of the canonical helium-burning locus, 
because the evolution accelerates towards the latter stages of 
the HB phase (i.e., from the bottom to the top of the canonical 
HB locus). Because we use the BHB as a standard candle to 
normalize the models, we can make an appropriate compari- 
son of the observed EHB locus and the predicted EHB locus. 

NGC 2419 We have assumed a metallicity of [Fe/H] = 
-2.12 (Harris 1996) and a distance of 83.2 kpc (Ripepi et 
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al. 2007). For the reddening, our assumption of E(B — V) = 

0. 03 mag is somewhat lower than that assumed elsewhere 
(0.08-0.1 1 mag; Harris 1996; Piotto et al. 2002; Ripepi et al. 
2007), but the upward shift in the models is required to align 
the observed and theoretical loci at the BHB. As one pro- 
gresses blueward in m/r[/v-»JA't/v color from l.Oto-0.8 mag, 
the observed HB appears to slope slightly upward relative to 
canonical HB, but it is worth noting that optical CMDs span- 
ning the B, V, and / bandpasses (Piotto et al. 2002; Ripepi et 
al. 2007) do not show a sloping HB that would justify nor- 
malizing the canonical locus significantly below the BHB. As 
one moves to stars bluer than the BHB, the observed HB lo- 
cus drops abruptly to fainter magnitudes, such that nearly all 
of the observed EHB stars are subluminous (i.e., BH stars). 

NGC 6273 We have assumed a metalUcity of [Fe/H] = 
-1.68 and a distance of 8.7 kpc (Piotto et al. 2002). For the 
reddening, our assumption of E{B-V) = 0.43 mag is larger 
than that assumed by Piotto et al. (2002; 0.37 mag), shift- 
ing the theoretical locus fainter by ^0.5 mag in the far-UV to 
give alignment with the observed BHB. As with NGC 2419, 
the HB slope appears somewhat steeper than expected as one 
moves blueward from 0.5 mag to -0.5 mag in mpuv—mi^uv 
color, even though the optical CMD (e.g., Piotto et al. 2002) 
shows no indication of an upward HB slope from the red HB 
to the BHB. At the hot end of the HB, the observed population 
is roughly split between normal EHB stars and BH stars. 

NGC 6715 We have assumed a metalUcity of [Fe/H] = 
-1.58 and a distance of 26.8 kpc (Harris 1996). For the red- 
dening, our assumption of E{B-V) = 0.16 mag is slightly 
higher than that assumed elsewhere (0.15 mag; Harris 1996; 
Siegel et al. 2007), shifting the models fainter by 0.08 mag to 
align with the observed BHB. The observed HB locus and the 
theoretical locus agree well over the entire BHB range, but 
at the EHB, the observed population is roughly split between 
normal EHB stars and BH stars. Note that NGC 6715 may 
actually represent the core of the disrupted Sgr dwarf galaxy 
(e.g., Layden & Sarajedini 2000; Siegel et al. 2007). Even if 
it is not the core of the Sgr dwarf, the earlier optical CMDs 
in this vicinity exhibit populations more complex than one 
would usually associate with a globular cluster, which may 
explain the large scatter in our UV CMD. 

NGC 2808 This cluster was examined extensively in Paper 

1. The cluster is included here for comparison, and we assume 
the same cluster parameters used previously. Although the 
HB sample in NGC 2808 is smaller than that in NGC 6715, 
the morphology of the observed HB locus is similar. The sim- 
ilarities between the NGC 6715 and NGC 2808 CMDs may 
be related to the similarities in metallicity and central con- 
centration; NGC 6273 is also close in metallicity, but is less 
centrally concentrated and less luminous than NGC 6715 and 
NGC 2808 (Harris 1996). 

NGC 6388 We have assumed a metalUcity of [Fe/H] = -0.6 
and a distance of 11.5 kpc (Piotto et al. 2002). For the red- 
dening, our assumption of E{B-V) = 0.41 mag is slightly 
larger than that of Piotto et al. (2002; 0.4 mag), placing the 
canonical locus 0.5 mag fainter than the observed BHB, as 
found by Busso et al. (2007) and discussed above. How- 
ever, Busso et al. (2007) assumed a distance of 11.2 kpc and 
E(B-V) = 0.45 mag; if we adopted these values, the canonical 
HB locus would be an additional 0.3 mag fainter than we have 
assumed. In contrast, assuming the shorter distance (10 kpc) 
and lower reddening (0.37 mag) of Harris (1996) would aUgn 
the canonical locus with the observed HB locus at the BHB, as 
we did for the metal-poor clusters. Although we have chosen 



an extinction that places the canonical BHB 0.5 mag below 
the observed BHB, nearly all of the observed EHB stars are 
subluminous with respect to the theoretical locus, as also seen 
in NGC 2419. However, unlike the situation in NGC 2419, 
the subluminous EHB stars of NGC 6388 cover a much wider 
range of color than that seen in any of the metal-poor clusters; 
the subluminous EHB stars in NGC 2419 are well-separated 
in color from the BHB stars, but there is no such color gap in 
the NGC 6388 population. Despite their extension to unusu- 
aUy red colors, the subluminous EHB stars in NGC 6388 are 
much bluer than the BSS (see Figure 3), and thus there is no 
confusion between these two classes of stars. Moreover, the 
subluminous EHB stars in NGC 6388 comprise a contiguous 
and distinct group, which suggests that they may aU be BH 
stars, including the reddest stars in this group. 

NGC 6441 We have assumed a metalUcity of [Fe/H] = 
—0.53, and a distance of 1 1.7 kpc (Harris 1996). For the red- 
dening, our assumption of E{B-V) = 0.53 mag is larger than 
that of Harris (1996; 0.47 mag), placing the canonical locus 
0.5 mag fainter than the observed BHB, as found by Busso 
et al. (2007) and discussed above. However, assuming the 
distance (15.2 kpc) and reddening (0.48 mag) of Busso et al. 
(2007) would place the canonical locus an additional 0.2 mag 
fainter than we have assumed. As with NGC 6388, assum- 
ing the Harris (1996) distance and reddening would align the 
canonical locus with the observed HB locus on the BHB, 
which is the normaUzation we used with the metal-poor clus- 
ters. Piotto et al. (2002) assumed a shorter distance (1 1 .3 kpc) 
and less reddening (0.44 mag); assuming these values would 
place our models above the observed HB locus for all tem- 
peratures (i.e., not simply the potential BH stars). Regardless 
of the normalization, the colors imply that there are almost 
no stars hotter than ~ 15,000 K in the cluster on or below the 
canonical EHB. There are stars below the canonical HB in the 
temperature range of 10,000 - 15,000 K, which is hotter than 
the BSS. When considered in the context of the NGC 6388 
BH population (which extends far to the red), it is possible 
that these subluminous stars in NGC 6441 are BH stars, de- 
spite their relatively red color. We inspected archival WFPC2 
images in the F336W (U) and F439W (B) bands to determine 
where these unusually red BH candidates would fall in optical 
CMDs. Although the core of the cluster is crowded at these 
longer wavelengths, five such stars were sufficiently isolated 
to obtain clean photometry, and in the B vs. U-B CMD they 
clearly fall within the low-luminosity end of the long blue HB 
tail, reinforcing the need for UV photometry to distinguish 
between temperature and luminosity variations in hot stars. 

In all six clusters, the hottest HB stars observed are sublu- 
minous with respect to canonical models (see Figures 3 and 
5). Thus there are BH stars in each cluster, with the possible 
exception of NGC 644 1 , where these subluminous stars are all 
relatively red. It is clear that canonical models at 7 = 0.23 can- 
not explain the HB populations of these clusters. Regardless 
of how such models are normalized via varying the assumed 
reddening or distance, they cannot reproduce the observed lu- 
minosity difference between the BHB and the subluminous 
EHB stars. 

4.2. Canonical Models with Enhanced Helium Abundance 

We next consider models with the cluster metallicity but 
with the helium abundance increased to F = 0.4. This value 
of helium enhancement is at the upper end of the plausible 
range for stars that have been enhanced at birth due to cluster 
self-enrichment. In Figure 6 we compare these models with 
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Y = 0.4 to those with Y = 0.23, assuming the same cluster 
parameters (|Fe/H], distance, and reddening) used previously. 

Figure 6 shows that an enhanced helium abundance pro- 
duces a more pronounced dropoff in luminosity between the 
BHB and the EHB. Relative to the Y = 0.23 locus, the F = 0.4 
locus is fainter on the EHB but much brighter on the BHB. 
The fainter EHB luminosity of the Y = 0.4 locus is due to a 
slightly smaller core mass while the brighter BHB luminosity 
is due to a greater energy output of the H-burning shell. This 
dependence of the luminosity difference between the BHB 
and EHB on the helium abundance is well established and 
was evident even in the early HB calculations of Sweigart & 
Gross (1976). We note in Figure 6 that the Y = 0.4 models 
for the very metal -poor cluster NGC 2419 are restricted to 
hot temperatures, even for those HB stars that lost no mass 
on the RGB. This is because the main sequence tumoff mass 
for a population at F = 0.4, [Fe/H]=-2.12, and an age of 13 
Gyr is only 0.614 Mq; this low tumoff mass, combined with 
this metal-poor composition, will only produce blue HB stars, 
even with no mass loss. 

It is clear that the F = 0.4 models do not agree well with the 
observed photometry when we assume the same parameters 
([Fe/H], distance, and reddening) used in the Y = 0.23 mod- 
els. The hot end of the F = 0.4 locus is not as faint nor as red 
as the observed BH stars, and the cooler part of the F = 0.4 
locus is much brighter than the observed BHB stars. While a 
combination of the F = 0.23 and F = 0.4 models would lead 
to a larger luminosity width for the EHB, the predicted range 
in luminosity would still not extend to the faint magnitudes 
reached by the BH stars. It is also clear that any range in F 
along the BHB must be fairly small, given the narrow lumi- 
nosity width of the BHB loci in Figure 6. Note that the larger 
scatter in NGC 6715 may be due to the complex populations 
associated with the Sgr dwarf galaxy (Lay den & Sarajedini 
2000; Siegeletal. 2007). 

At first glance, one might think that the agreement between 
F = 0.4 locus and the observed BHB and EHB stars could 
be improved by assuming a larger reddening or distance, be- 
cause such adjustments to the cluster parameters would shift 
the models to fainter magnitudes. Indeed, if the Y = 0.4 lo- 
cus were shifted until it fit the observed BHB, the lower edge 
of the EHB for these models would be in much better agree- 
ment with the observed lower edge of the BH stars. However, 
such adjustments to the reddening or distance would violate 
important observational constraints. 

If we attempt to fit the BHB with the F = 0.4 locus, we wiU 
have to confront the following observational difficulties. Opti- 
cal CMDs of NGC 2419 (Sandquist & Hess 2008), NGC 6273 
(Piotto et al. 1999), NGC 6715 (Layden & Sarajedini 2000), 
or NGC 2808 (Bedin et al. 2000) do not show a significant 
difference in the V-band luminosity of the red HB and BHB 
stars, implying that any difference in the helium abundance 
between these stars must be small. Thus the cluster param- 
eters we derived in fitting the BHBs of these clusters to the 
canonical F = 0.23 locus should also apply to the red HB 
stars. In the case of the two metal-rich clusters NGC 6388 
and NGC 6441, it is well known that the HB slopes upward 
(Rich et al. 1997; Busso et al. 2007), with the V-band lumi- 
nosity of the BHB being ~0.5 mag brighter than that of the 
red HB clump. The BHB stars in these clusters may thus have 
an intermediate helium abundance (F ~ 0.3). However, as dis- 
cussed above, we took this offset in the BHB luminosity into 
account when fitting the F = 0.23 model to our UV CMDs. 
Thus our parameters for these two clusters should also apply 



to their red HB clumps. If the BHB should be dominated by 
stars with F = 0.4 and the red HB by stars with F = 0.23, one 
would therefore predict a luminosity difference between the 
BHB and red HB clump that is much larger than observed 
in any of our six clusters. To avoid this contradiction, one 
might argue alternatively that the red HB stars are themselves 
helium-rich, but this would imply a much brighter red HB lu- 
minosity and hence a much younger (and untenable) cluster 
age from the observed luminosity difference between the red 
HB and the MSTO (e.g., see Martin-Franch et al. 2010). As 
one example, Catelan et al. (2006) have compared the CMDs 
of NGC 6388 and 47 Tuc and concluded that both of these 
clusters have similarly old ages. 

In short, while the F = 0.4 locus exhibits a sharp drop in 
luminosity between the BHB and EHB, these models cannot 
by themselves explain the BH phenomenon. Given the lack 
of helium enhancement on the red HB and the relative lumi- 
nosities of the BHB and red HB, we cannot align models of 
arbitrarily high F to the BHB without violating other observa- 
tional constraints. If we assume the BHB has a standard he- 
lium abundance (F = 0.23, in the case of the four metal-poor 
clusters) or an intermediate helium abundance (F ~ 0.3, in the 
case of the two metal-rich clusters), then even a wide range of 
helium abundance (0.23 < F < 0.4) on the EHB would not 
explain either the faint luminosities or the wide color range 
observed in the BH stars. 

4.3. Flash Mixed Models 

We next consider models which undergo flash mixing due 
to high mass loss on the RGB. In Figure 7, we plot these flash- 
mixed models in our UV CMDs, assuming the same redden- 
ing and distance values obtained in section 4. 1 . Those values 
align a canonical F = 0.23 locus to the BHB in the four metal- 
poor clusters, but place this Y = 0.23 locus 0.5 mag fainter (in 
nifuy) than the BHB in the two metal-rich clusters. For each 
cluster in Figure 7 we show the HB evolutionary tracks for 
three representative flash-mixed models; two of these tracks 
faU at the extremes of the range in the mass loss rate leading 
to flash mixing, as measured by the value of while one 
falls halfway between these extremes. We plot these triplets 
of flash-mixed tracks for both F = 0.23 and F = 0.4 in or- 
der to bracket the range of behavior expected from any sub- 
population in the cluster that undergoes enough mass loss for 
flash mixing. All of the flash-mixed tracks for each cluster 
fall within a narrow temperature range. The offset in lumi- 
nosity among these tracks is due to small differences in the 
helium-core mass. 

As can be seen in Figure 7, flash-mixed models can repro- 
duce the faint luminosities observed in the BH population, 
particularly the flash-mixed models originating in a subpopu- 
lation with enhanced helium. Flash mixing thus offers a way 
of reproducing the faint luminosities observed in the BH stars 
without causing problems at the BHB or red HB clump, as dis- 
cussed in the last section. Furthermore, the flash-mixed mod- 
els, when combined with the canonical HB locus, are better 
able to reproduce the liirge luminosity width observed at the 
hotendof flieEHB. 

There is one obvious problem, however, which was also 
seen in the previous sections; like canonical models, flash- 
mixed models (whether originating in a population with en- 
hanced or standard F) cannot reproduce the red colors in the 
BH population. Paper I explored various ways of producing 
such red BH stars, including enhanced line blanketing in the 
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atmosphere, due to an increase in the atmospheric Fe abun- 
dance via radiative diffusion. In the case of NGC 2808, an 
increase in the [Fe/H] from -1.36 to +1.0 could shift the BH 
stars -^0.3 mag to the red. However, this color shift is not 
enough to explain the full color range of the BH population. 
Furthermore, the effect would be even smaller for the metal- 
rich clusters, which are the very clusters exhibiting the red- 
dest BH stars. Because the metal-rich clusters already incur 
significant line blanketing, increasing the [Fe/H] via diffusion 
does not help much; for example, increasing the [Fe/H] in 
NGC 6388 from -0.6 to -l-l .4 only shifts the BH color redward 
by 0.26 mag. Thus, at this time, the mechanism for producing 
subluminous KB stars with such red colors remains a puzzle, 
and it suggests that the formation mechanism for these stars 
is not yet fully understood. 

5. DISCUSSION 

5.1. Blue Hook Stars 

In this paper, we have compared the UV CMDs of six 
massive globular clusters with the expectations from stellar 
evolution theory. We find that in each of these clusters, the 
hottest HB stars are subluminous with respect to canonical 
models, whether or not those models assume a standard he- 
lium abundance (Y = 0.23) or an enhanced helium abundance 
(F = 0.4). The most straightforward interpretation is that all 
of these clusters host BH stars. However, in the two metal- 
rich clusters, the subluminous stars either extend to relatively 
red colors (NGC 6388) or are restricted to these red colors 
(NGC 6441), which might imply that an additional or dis- 
tinct phenomenon is at work at high metallicity. Indeed, 
the existence of BH stars in NGC 6441 has been questioned 
(e.g., Dieball et al. 2009). The BH population of NGC 6388 
may be the link between the BH populations in the metal- 
poor clusters and that in NGC 6441, because it spans the full 
color range observed in these clusters. The bluest stars in the 
NGC 6388 BH population appear to be analogs of the BH 
stars in the metal-poor clusters, while the reddest stars in the 
NGC 6388 BH population appear to be analogs of the BH 
stars in the NGC 6441 population. Moreover, the entire BH 
population in NGC 6388 appears to be contiguous. 

In some massive clusters, all of the EHB stars may have 
originated in a sub-population with enhanced heUum. For ex- 
ample, in NGC 2808, nearly all of the EHB stars falling within 
the Y = 0.23 canonical locus (Figure 6; pink shaded region) 
also fall within the 7 = 0.4 canonical locus (Figure 6; green 
shaded region); such stars are normal in the sense that they 
are not flash-mixed. However, there is a dearth of stars falUng 
in the brighter part of the Y = 0.23 EHB locus, which may 
indicate that the dominant Y = 0.23 population in this cluster 
only produced BHB and red HB stars. In NGC 2419, the sit- 
uation is even more extreme. Few EHB stars fall in either the 
Y = 0.23 or y = 0.4 canonical loci, and thus all of the EHB 
stars in this cluster may have been produced via flash-mixing 
(see Figure 7). As with NGC 2808, the dominant Y = 0.23 
population in NGC 2419 may only have produced BHB and 
red HB stars. 

5.2. The Case for Flash Mixing 

The faint luminosities observed in the BH stars can be re- 
produced by models that undergo flash mixing during a late 
heUum-core flash on the WD cooling curve. Helium enhance- 
ments may also play an important role in the formation of BH 
stars, primarily because a higher helium abundance lowers the 



MSTO mass at a given age so that less mass loss is required on 
the RGB to populate the hot HB . Helium enhancements also 
depress the luminosity somewhat in canonical EHB models, 
but not by the ~1 mag of the faintest BH stars. 

Several recent studies have argued against flash mixing as 
the origin of the BH stars. For example, Lee et al. (2005) 
attribute the BH stars to the evolution of a population with 
a high initial helium abundance (Y = 0.4). However, as we 
demonstrate in §4.2, an enhancement in Y cannot by itself ex- 
plain the BH phenomenon. While an enhanced Y decreases 
the EHB luminosity and increases the BHB luminosity (thus 
providing a hook at the hot end of the HB), one cannot arbi- 
trarily adjust the cluster reddening and distance to match the 
theoretical HB locus to an observed CMD without regard to 
the relation between the BHB, the red HB, and the MSTO. 
Furthermore, the models of Lee et al. (2005) do not predict 
any carbon enhancement in the BH stars, while optical spec- 
tra of BH stars in u) Cen and NGC 2808 do suggest such an 
enhancement may be present (Moehler et al. 2004, 2007) - an 
enhancement naturally produced by flash-mixing. 

D' Antona et al. (2010) go further, postulating that MS stars 
born at F « 0.4 could undergo deep mixing on the RGB that 
elevates their helium abundance to 7 w 0.8 before arrival on 
the EHB. The mechanism that drives this deep mixing is un- 
clear, unlike the mechanism for flash mixing that follows from 
the basic stellar structure equations and which is the inevitable 
outcome of a helium-core flash on the WD cooUng curve. 
However, the effect of the D' Antona et al. (2010) "RGB mix- 
ing" is to generate EHB stars with very high helium abun- 
dances that are hotter than the canonical HB, and in this way 
their models are similar to those produced by flash mixing. 
An important distinction, however, is that their RGB mixing 
precedes the helium-core flash, and thus does not produce an 
enhancement in the carbon abundance. 

Castellani et al. (2006) have argued that canonical HB mod- 
els can reach the high Teft observed for BH stars, but their 
EHB stars reach a maximum Teff of ~32,000 K. This is much 
less than that found by Moehler et al. (2004, 2007) for the 
BH stars (>35,000 K), and much less than the temperatures 
reached by either flash mixing or RGB mixing models. Fur- 
thermore, canonical HB models would not produce any en- 
hancements in helium or carbon, as observed for such stars 
(Moehler et al. 2004, 2007). 

In Paper 11, we found that a delayed helium-core flash could 
sometimes lead to a more shallow form of flash mixing, where 
the envelope hydrogen is only mixed into the outer layers of 
the core. When this happens, flash-mixed stars can retain 
some residual envelope hydrogen, and consequently would 
be distributed closer to the canonical EHB. Our models im- 
ply that shallow mixing may be a relatively rare occurrence, 
especiaUy at low metalUcities. Our CMDs do not show any 
systematic variation in the number of BH stars falling close 
to the canonical locus, suggesting that shallow mixing is not 
important in the metallicity range covered by our clusters. It 
may be, however, that we simply do not have enough BH stars 
to detect such a trend. With the current sample, it appears 
that the BH stars in each cluster fall up to ~1 mag below the 
canonical HB, regardless of metaUicity. 

The sdB and sdO stars of the Galactic field are the analogs 
of the EHB and post-EHB stars found in clusters. One would 
expect some of these stars to also show evidence of flash 
mixing. This expectation motivated far-UV spectroscopy of 
heUum-rich sdB stars in the Galactic field using the Far Ultra- 
violet Spectroscopic Explorer, to see if these stars exhibit the 
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helium and carbon abundance enhancements expected from 
flash mixing. Those observations do in fact lend support 
to the flash mixing scenario. The spectrum of PG1544+488 
(archetype of the He-sdB class) shows incredibly strong car- 
bon lines, indicating a surface that is » 2% carbon by mass 
(Paper II). 

Further insight into the nature of the BH stars may come 
from upcoming far-UV spectroscopy of NGC 2808. In 
the current observing cycle, we will use the Space Tele- 
scope Imaging Spectrograph on HST to obtain spectra of BH 
stars, normal EHB stars, and BHB stars selected from the 
NGC 2808 CMD (Figure 5). With these spectra, we will be 
able to search for the atmospheric signatures of flash mixing 
(e.g., enhanced carbon, diminished H) and perhaps shed light 
on the origin of the wide color range exhibited in BH stars. 

Our understanding of the BH phenomenon is far from com- 
plete. None of the models discussed in this paper (canonical at 
standard 7, canonical at enhanced 7, or flash-mixed) are able 
to explain the red colors found in some or all of the sublumi- 
nous stars in each cluster. This is particularly true in the two 
metal-rich clusters in our sample. Furthermore, the near-total 
lack of normal EHB stars in half of our clusters is puzzling. 
If we consider the categories of HB stars (red HB, BHB, nor- 
mal EHB, BH) as a sequence governed by cluster parame- 
ters (mass loss, helium abundance, cluster mass, etc.), why do 
some clusters (NGC 2419, NGC 6388, and NGC 6441) ex- 
hibit BHB stars and BH stars but few normal EHB stars? The 
dearth of such stars may imply a gap in one or more cluster 
parameters, but the mechanism for such a gap is unclear. 

5.3. Why Massive Clusters? 

A remarkable feature of all clusters either known or sus- 
pected to host BH stars is that they are among the most mas- 
sive clusters in the Galaxy. DiebaU et al. (2009) also noted 
that BH stars only appear to exist in massive clusters, and sug- 
gested that this could simply be a statistical effect, because 
"intrinsically rare objects can only be found by searching a 
sufficiently large number of stars." To examine this possibil- 
ity, DiebaU et al. (2009) apply a statistical test by comparing 
the locations of clusters with and without significant BH pop- 
ulations within the cumulative mass distribution of the clus- 
ters in their sample. Although Dieball et al. (2009) state that 
their analysis is "effectively combining many low-mass clus- 
ters to form aggregate high-mass clusters," their test, in fact, 
ignores the actual numbers of BH stars in each cluster, and 
instead only considers whether clusters do or do not host sig- 
nificant numbers of BH stars. This is a difficult test, given 
the small number of clusters in their comparison. Although 
Dieball et al. (2009) suggested that it would be worthwhile 
to investigate how the numbers of BH stars scale with cluster 
properties, they did not actually perform this test due to the 
non-uniformities in the extant data. A uniform data sample 
would be important if one were trying to address the total size 
of the BH population as a function of cluster mass. However, 
the Dieball et al. (2009) sample already contains sufficient 
data to explore the relationship between BH population size 
and cluster mass. If the UV imaging is sufficient to produce 
a statistically significant sample of EHB stars in each cluster, 
one can ask if a significant fraction of the EHB stars are BH 
stars, and can then compare this fraction in clusters of differ- 
ent masses. This comparison is meaningful even if different 
clusters are observed in distinct ways (e.g., differences in re- 
gion sampled, bandpasses used, depth); the important point is 
to observe a large enough sample of EHB stars to determine 



the relative numbers of normal and subluminous (BH) EHB 
stars. 

We refer to Table 2 of DiebaU et al. (2009). Their ta- 
ble Usts 5 clusters hosting significant numbers of BH stars 
(NGC 2419, NGC 2808, w Cen, NGC 6388, and NGC 6715). 
All of these clusters have masses ranging from 1.42 x 10'' Mq 
(NGC 2808) to 3.35x10*' Mq {lo Cen). We would add 
NGC 6273 to this Ust; because it is 0.2 mag fainter in My 
than NGC 2808 (Harris 1996), its mass is «1.2x 10^ Mq. Ta- 
ble 2 of Dieball et al. (2009) also lists 11 clusters that do not 
have significant numbers of BH stars. Of course, besides these 
1 1 clusters, there are many more low-mass and intermediate- 
mass clusters without known BH stars, but these 1 1 were con- 
sidered by Dieball et al. (2009) to have sufficient imagery to 
rule one way or the other. Among that list of 1 1 clusters is 
NGC 6441, with a mass of 1.57x 10'' Mq; we would move it 
to the list of clusters hosting BH stars. Also in that list of 1 1 
clusters is NGC 104 (47 Tuc); it is massive (1.5x10'' Mq), 
but it is the only cluster in the entire table that has no BHB 
stars (its HB population is entirely red). The total mass of 
the 9 remaining clusters in this list is 5.1 x lO*" Mq. The stel- 
lar mass in this ensemble of 9 clusters exceeds the mass in 
any individual cluster, and is twice the mass of any cluster 
we have explored in our own study. Although this ensem- 
ble is comprised of 9 clusters hosting BHB stars and in some 
cases (e.g., M13, M80, NGC 6752) substantial populations of 
EHB stars, this ensemble has few if any BH stars (possibly 
3-5; DiebaU et al. 2009; Sandquist et al. 2010). If finding BH 
stars simply required the searching of enough stellar mass, we 
would expect to find tens of BH stars in the ensemble mass of 
9 clusters totaling 5.1x10'' Mq. The impUcation is that any 
cluster hosting BHB stars will also host BH stars if its mass 
exceeds ~1.2x lO'' M©. However, the fraction of EHB stars 
that are subluminous varies from cluster to cluster in our sam- 
ple, and that fraction does not appear to correlate with either 
cluster mass or metallicity, implying additional factors are at 
work (including perhaps the dynamical history of interactions 
between the cluster and the Galaxy). In half of our sample 
(NGC 2419, NGC 6388, and possibly NGC 6441), virtually 
all of the EHB stars appear to be BH stars, with very few 
stars on the canonical EHB, whereas in the remaining three 
clusters (NGC 6273, NGC 6715, and NGC 2808), the spUt 
between normal and subluminous EHB stars appears to be 
roughly equal. 

Note that low-mass clusters can have large samples of EHB 
stars without having any BH stars. Consider NGC 6752, 
which is one of the 9 clusters discussed above that hosts no 
BH stars; it is Usted in Table 2 of DiebaU et al. (2009) witii 
a mass of 0.32x10'' Mq. For NGC 6752 (Landsman et al. 
1996) and NGC 2808 (Paper I), there exist CMDs that in- 
clude a far-UV bandpass. In the case of NGC 6752, the CMD 
is sampling a less massive cluster outside of the core, while in 
the case of NGC 2808, the CMD is sampUng a more massive 
cluster within the core. Despite these distinctions, the samples 
of EHB stars and their descendents are comparable in size. 
In each CMD, there are --80 HB stars hotter than 16000 K, 
along with fourpost-EHB stars - i.e., the shorter-lived AGBM 
stars. However, for the case of NGC 2808, 46 of the EHB stars 
are BH stars (lying weU below the canonical EHB), while in 
NGC 6752, none of them are. Given that the size of the EHB 
samples are approximately the same, the number of BH stars 
in the two clusters are discrepant at a high statistical signif- 
icance. This argument is not restricted to NGC 6752; other 
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low-mass clusters, such as M13 and M80 (Ferraro et al. 1998) 
host dozens of EHB stars but apparently few if any BH stars 
(M13 may host 2 BH stars; Sandquist et al. 2010). The large 
numbers of EHB stars in both the low-mass cluster sample 
and the high-mass cluster sample provide strong support for 
the notion that the fraction of EHB stars that are subluminous 
(i.e., BH stars) is nearly zero in low-mass clusters and ^50- 
100% in high-mass clusters. 

There are several possible reasons why BH stars may only 
appear in massive globular clusters. One possibility is that 
only massive clusters can retain the helium-rich ejecta from 
their first generation of stars. Besides this helium enrichment, 
such self-pollution might account for many of the abundance 
variations in O, Na, Mg and Al observed in MS and RGB stars 
(see Gratton et al. 2004 and references therein). 
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An alternative hypothesis is that massive clusters have 
many close binaries, which produce the high mass loss needed 
for flash mixing. Of the massive clusters known to host BH 
stars, none are core-collapsed, and all are of intermediate cen- 
tral concentration (Harris 1996), so that one might expect a 
significant fraction of surviving binaries (Ivanova et al. 2005). 
However, the lack of a radial gradient in the NGC 2808 BH 
population argues against a binary origin (Bedin et al. 2000). 



Support for proposal 10815 is provided by NASA through 
a grant from STScl, which is operated by AURA, Inc., un- 
der NASA contract NAS 5-26555. We thank the anonymous 
referee for comments that helped improve the clarity of this 
paper. 



Harris, W.E. 1996^ AJ, 112, 1487 

Heyer, I., Richardson, M., Whitmore, B., Lubin, L. 2004, WFPC2 Instrument 

Science Report WFPC2 2004-001 
Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875 

Ivanova, N., Belczynski, K., Fregeau, J.M., & Rasio, FA. 2005, MNRAS, 
358, 57 2 

Landsman, W.B., Sweigart, A.V, Bohlin, R.C, Neff, S.G., O'Connell, R.W., 
Roberts, M.S., Smith, A.M., & Stecher, T.P 1996, ApJ, 472, L93 

Lanz, T., Brown, T.M., Sweigart, A.V, Hubeny, I., & Landsman, W.B. 2004, 
ApJ, 602, 342 (Paper II) 

Lanz, T., & Hubeny, 1. 2003, ApJS, 146, 417 

Lanz, T., & Hubeny, 1. 2007, ApJS, 169, 83 

Layden, A.C, & Sarajedini, A. 2000, AJ, 119, 1760 

Lee, Y.-W., Joo, S.-J., Han, A.-L, Chung, C, Ree, C.H., Sohn, Y.-J., Kim, 

Y.-C, Yoon, S.-J., Yi, S.K., & Demarque, P 2005, ApJ, 621, L57 
Martin-Franch, A., Cassisi, S., Aparicio, A., & Pietrinfemi, A. 2010, ApJ, in 

press, astro-ph/1003.5489 
Maybhate, A., et al. 2010, 'ACS Instrument Handbook," Version 9.0 

(Baltimore: STScI). 
Miller Bertolami, M.M., Althaus, E.G., Unglaub, K., & Weiss, A. 2009, 

JPhCS, 172, 12014 
Milone, A.P, et al. 2008, ApJ, 673, 241 

Moehler, S., Dreizler, S., Lanz, T., Bono, G., Sweigart, A.V, Calamida, A., 

Monelli, M., & Nonino, M. 2007, A&A, 475, L5 
Moehler, S., Sweigart, A.V, Landsman, W.B., Hammer, N.J., & Dreizler, S. 

2004, A&A, 415, 313 
Piotto, G., Bedin, L.R., Anderson, J., King, I.R., Cassisi, S., Milone, A.P, 

Villanova, S., Pietrinfemi, A., & Renzini, A. 2007, ApJ, 661, L53 
Piotto, G., et al. 2002, A&A, 391, 945 
Piotto, G., et al. 2005, ApJ, 621, 777 

Piotto, G., Zoccali, M., King, I.R., Djorgovski, S.G., Sosin, C, Rich, R.M., 

& Meylan, G. 1999, ApJ, 118, 1727 
Reimers, D. 1975, Mem. Roy. Soc. Liege, 8, 369 
Reimers, D. 1977, A&A, 57, 395 

Renzini, A. 1990, in ASP Conf. Ser. 11, Confrontation Between Stellar 
Pulsation and Evolution, ed. C. Cacciari & G. Clementini (San Franciso: 
ASP), 549 

Renzini, A. 2008, MNRAS, 391, 354 

Rich, R.M., et al. 1997, ApJ, 474, L25 

Ripepi, v., et al. 2007, ApJ, 667, L61 

Rosenberg, A., Recio-Blanco, A., Garca-Marii, M. 2004 ApJ, 603, 135 

Salaris, M., Chieffi, A., & Straniero, O. 1993, ApJ, 414, 580 

Sandquist, E.L., Gordon, M., Levine, D., & Bolte, M. 2010, AJ, 139, 2374 

Sandquist, E.L., & Hess, J.M. 2008, AJ, 136, 2259 

Sirianni, M., et al. 2005, PASP 117, 1049 

Siegel, M.H., et al. 2007, ApJ, 667, L57 

Sweigart, A.V. 1997, in The Third Conference on Faint Blue Stars, ed. A.G.D. 

Phillip, J. Liebert, & R.A. Saffer (Schenectady: L. Davis Press), 3 
Sweigart, A.V, & Gross, PG. 1976, ApJS, 32, 367 



Updated at |http://physwww.physics.mcmaster.ca/~harris/inwgc.dat 



